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The disturbances produced by a slightly inhomogeneous random medium on a passing 
wave can be classified into contributions depending on an increasing number of successive 
scatterings. The individual contributions appear in an expansion of the solution of an in- 
tegral equation. The first term, the Born approximation, only accounts for a single scat- 
tering. Convenient expressions for this approximation result from a saddle point treatment for 
short distances, and from a Fraunhofer approximation for larger distances. The evaluation 
of the higher-order contributions, describing plural-scattering effects, leads to mathemat- 
ical difficulties which are evaded by considering the scattering mechanism as a Markovian 
process. The corresponding theory can be developed with the aid of an integro-differen- 
tial diffusion equation; the latter refers to the joint probability density of the lateral and 
angular deviations suffered by the trajectory of the passing wave. The equation in question 
can be solved with the aid of four-dimensional operational calculus; it reduces to the simple 
differential equation of Fokker-Planck under special conditions. The application of the gen- 
eral equation to tropospheric point-to-point radio communication is worked out. It is shown 
that the far-distance field, associated in this case with multiple scattering, does decrease 
proportionally to the second or third power of the inverse distance. 

1. Introduction 

By volume scattering we understand the property that each volume-element of a slightly 
inhomogeneous medium does deviate a fraction of the energy of any incident radiation from its 
original propagation direction. In its simplest form such scattering is described by a scalar 
wave equation with a term consisting of the product of the wave function and the fluctuating 
component of the refractive index. By solving this equation while treating this term as if it 
were known w^e obtain an integral equation. The successive terms of the Neumann expansion 
of its solution represent the primary field (corresponding to a medium without a random com- 
ponent of the refractive index), the contribution due to a single scattering, that produced by 
two successive scatterings, three successive scatterings, and so on. 

Most attention has been paid to the single-scattering contribution which is known as the 
Born approximation. Convenient expressions can be deduced from its rigorous representation, 
either by applying a saddlepoint method, or by taking account of the dimensions of the relevant 
coherently scattering parts of space ("blobs") which are small compared to their distances to 
the transmitter and the receiver. The saddlepoint method involves geometric-optical approxi- 
mations for the distortions of the wave fronts which would leave the transmitter undisturbed 
in a nonfluctuating medium. The other method (Fraunhofer approximation) leads to the 
distribution of the scattered energy over the various directions around the original propagation 
direction. In either method the correlation between the field strengths observed at two 
different (or identical) places at two different (or identical) moments can be derived from a 
corresponding correlation function for the random fluctuations of the refractive index. 

Mathematical difficulties till now prevented a similar treatment of the higher-order 
scattering contributions. Therefore, the successive scatterings have simply been assumed as 
completely independent of each other, which is certainly justified if they do not take place too 
close to each other. This Markovian scattering process can be described with the aid of a joint 
probability density h z which combines the chances for special, both lateral and angular, devia- 
tions away from the rectilinear propagation path covered in the case of a homogeneous medium. 
This probability density satisfies a transport or diffusion equation which can be solved rigorously 
if the equation is first simplified in view of the assumption of small final angular deviations. 
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The equation in question reduces to a second-order Fokker- Planck equation if terms depending 
on the average fourth, sixth, etc., power of the angular deviations (connected with a single 
scattering) may be neglected. 

The probability density h z being known, it is possible to compute the field strength for 
tropospheric radio propagation if the radiation diagrams of both the transmitter and the 
receiver are given. The explicit evaluation is difficult but approximations for large distances, 
which include the effects of scatterings of any order, can be worked out for special models, for 
instance for Norton's modified Bessel model. Such calculations show the insufficiency of the 
corresponding expressions derived from the above-mentioned Fokker Planck equation. 

The various items indicated here are discussed in more detail in the subsequent sections. 

2. Integral Equation of a Slightly Inhomogeneous Medium 

The elimination of the magnetic field from Maxwell's equations for a medium with a die- 
lectric constant 

e=e + 5e(x,y, z) 

leads for time harmonic solutions proportional to e~ io}t to the equation: 

V 2 #+V(V In e • E) +a>WE=0. 

The gradient term can be neglected provided that the changes of 8e/e are small over one wave- 
length. The equation can then be put in the following form, for each component of E: 

V 2 E-\-oo 2 iJL eoE= — o) 2 iJL 8eE. 

By treating this wave equation as if the right-hand side were a known function, we get the 
following "solution" 

E(P)=EAP)+^j drMQ)E(Q) QP > (i) 

which, as a matter of fact, constitutes a scalar integral equation for each component of the 
unknown wave function E{P). The integration extends over all volume-elements dr Q of 
the space containing 8e fluctuations. The term E pr (P) represents the primary field, that is 
the solution in the absence of the random component 5e; QP is the distance from the integra- 
tion point Q to the point of observation P (receiving antenna) . 

The Neumann Liouville expansion of the solution of (1) starts with 

E(P)=EAP)+^rj dr Q 8e(Q)E pr (Q) -Qp 

Q+QP) 



+(l?) \ dT MQ)jdT Q MQ')E 9 rm ! QrQr l 

=E„(P)+E™(P)+E<»(P)+ . . ., say. (* =W/W. (2) 



The second term represents the contribution due to single scatterings and constitutes the Born 
approximation for the total scattered field. The third term E {2) (P) is recognized as the effect 
of two successive scatterings at two volume-elements dr Q and dr Q ^ ; the fourth term likewise 
represents the contribution of third-order scatterings, and so on. Most literature only con- 
cerns the Born approximation. 

3. Statistical Properties of the Born Approximation 

The random character of the fluctuations in slightly inhomogeneous media suggests to 
introduce statistical quantities right at the beginning, thus applying methods labelled as 
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"dishonest" in Keller's paper [Keller, 1962] on the subject. The statistical properties concern, 
e.g., the average field strength and the correlation between the fields observed at two neighbor- 
ing points Pi and P 2 , possibly at different moments U and t 2 . These properties can be derived 
from the quantity 

{E^{P u h)E^'(JP 2 ,t 2 )) 

for which we obtain, in view of the above expression for E {1 \ 

/oj 2 /xA 2 C pikoQiPi (* k f -ihQ2P 2 

(l?) J dT °L E »^ -QJ\) dr Q2 EUQ2) e -Q^(8e(Qi, tOSeiQ,, «). (3) 

A further evaluation needs an explicit statistical assumption concerning the average of 
the product of the dielectric-constant fluctuations at two different places. It is customary 
to introduce here the hypothesis of l 'homogeneous turbulence" according to which this average 
only depends on the relative positions of Q x and Q 2 , and the time difference t\— 1 2 . Neglecting 
further the time variations this amounts analytically to a dependence on the differences of the 
coordinates of Q x and Q 2 , or on the vector Q\Q 2 . The homogeneous turbulence can then be 
fixed completely with the aid of the following normalized autocorrelation function: 

CUQ1Q2) (8e2{Q)) • (4) 

We next assume a primary field due to a point source at T, and normalized according to the 
formula 

p ik TQ 
E »r(Q)=-fQ- (5) 

The substitution of (4) and (5) into (3) yields: 

/VunY f C e ifc (TQ l +Q l P 1 -rQ 2 -Q 2 P 2 ) > 

<^CP0*«'CPO>-(^) <^>J dr Ql j dr Qi TQi . QiPi . TQ2 . Q2 p 2 CUQM. (6) 

This expression with a double integration over the entire randomly fluctuating medium is basic 
for all investigations starting from the Born approximation. 

4. A Saddlepoint Method Applied to the Born Approximation 

Saddlepoint approximations are always applicable for frequencies which are sufficiently 
high. Its application to the Born approximation E (l) (P) as defined by (2) amounts to replacing 
the integration point Q in the denominator by its projection Q' on the line TP connecting the 
point source with the point of observation, and to expanding the exponential up to second- 
order terms with respect to the coordinates y q and z g ; the y and 2-axes are here assumed as 
perpendicular to the line TP constituting the z-axis. We then obtain, also using (5), 

B» CP) _-ft e«.-J<M«(<2) Q , T Q!p 

We pass from the dielectric constant e G -\-de to the corresponding expression 

n 2 =e/6 =l+28n 

for the refractive index n, so as to have de=2e 8n. A substitution of the Taylor expansion 
(up to second-order terms) for 8n Q as a function of Y Q and Z Q then leads to a result derived 
[Bremmer, 1958], which can be interpreted by a phase shift 

50=£ o ( Xp <fa Q '5n(Q'), (7) 
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and a relative change of the amplitude given by: 



8A 
A 



. 1 (Vr Q'T-Q'P (VH Vbn\ . 

" 2 J* r ^ TP W" + d«V e . (8) 



The phase correction is quite obvious in view of the fluctuations of the refractive index 
n. The amplitude change points to a lens effect. In fact, the regions with a negative or posi- 
tive value of (d 2 /()y 2 +c) 2 /dz 2 )5n cause a small increase or decrease of the field which is equiva- 
lent to a focusing or a defocusing. 

It is also possible to derive statistical averages for these geometric-optical quantities 
when assuming homogeneous turbulence. For instance, the variance (<5</> 2 ) of the phase 
fluctuations at a special point, or the correlation between these fluctuations at two different 
points are obtained from a double integral consisting of the product of two expressions (7), 
substituting the autocorrelation function corresponding to (4) for the refractive-index fluctua- 
tions. This autocorrelation function may be represented by 



Cbn (Ql Q2) = Csni j ) 



in which we introduce the scale of turbulence I in order to get a function of a dimensionless 
parameter. The variance in question then becomes as follows at a point at a distance TP=d 
from the point source, provided that d )) I: 

{bcf> 2 )=2kl{bn 2 )ld Jj dyCUy) =8d/d , say. (9) 

A similar computation, starting from (8) and to published elsewhere, leads to the following 
corresponding expression for the variance of the amplitude fluctuations: 

<(S)VMW-^}*<< 

Variances of other geometric-optical quantities, such as angular deviations from ray tra- 
jectories which are rectilinear in the absence of the random fluctuations, can be computed in a 
similar way. We refer in particular to the investigation [Muchmore and Wheelon, 1955] 
and to a discussion of the wave-front disturbances caused by turbulent random fluctuations 
[Bremmer, 1963]. 

5. Fraunhofer Approximations of the Born Approximation 

The limited range of noticeable coherence of the ^-fluctuations involves a splitting up of 
the integral of (6) into terms comprising a double integration over individual incoherently 
scattering volume elements. These elements or "blobs" may have dimensions of the order 
of the parameter I introduced above in the autocorrelation function Cdn- 

In the application to far-distance tropospheric propagation the effective blobs are situ- 
ated in a rather small part of space situated simultaneously above the tangenting plane to the 
earth (horizon plane) through the transmitter, and above that through the receiver. The re- 
sulting conventional theory is summarized below, emphasizing the influence of the finiteness 
of the blob dimensions. 

The size of the blobs being small compared to the distances Q\T, QzT, Q\P\, Q2P2 from their inner points to 
the transmitter and the points of observation, it is justified to replace these inner points by a fixed central 
point Qq of the blob in question when considering the denominator in (6) . Moreover, these distances may be 
approximated in the exponential by expressions which are linear with respect to the coordinates of Qi and #2 

(the origin being at Qo), that is the components of two vectors X\ and x 2 . We further introduce new mixed 

coordinates (components of two vectors £ and rj) according to: 
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£=a? 2 +£i, r} = x 2 —x 1 . 

In the case of spherical blobs of radius L/2 around their central point Q the splitting of (6) into contribu- 
tions over the various blobs then can finally be represented by: 



< g a,(P l)g a, ( p 2))= ((f^> g t ° qi Qapi Qop ^ J ^H,. Cl ,g). (11) 

Here we have introduced two other new vectors V and IT defined by: 



W = COS 2 U b — UTQ G , 

where Uq g p 2 marks a unit vector in the direction of Q0P2, the other unit vectors being defined likewise; b repre- 
sents the bissectrix of the angle \p = PiQ P 2 . 

The special case P { = P 2 leads at once to the differential scattering cross section a which represents the 
energy scattered by a unit scattering volume into a unit infinitesimal solid angle along the direction QqP, for 
a unit energy-current density of the incident radiation. This definition proves to be equivalent with that used 
later on (see sec. 7). The former definition implies the following value for the total scattering cross section 
connected with the N incoherently scattering blobs contained in a unit volume-element: 

tr=N TQl QoPHE^(P) E^*(P)) i; 

the index (1) here refers to the contribution connected with a single blob. Hence, in view of (11), 



JV(q,W 2 (fiC 2 ) 

128tt 2 

\i=H\<L 



J 4j4e<*° ( "««-"^Ca« (~y 



in which the new unit vectors are directed along the scattered and the incident radiation (along Q P and TQ {) ), 
respectively. 

The integral over £ amounts to the volume of the common space of two spheres with radius L and centers 
separated by a distance 2\y\. This common space only exists if \ri\<^L. The blobs (the spherical form of 
which has only been assumed in order to show the influence of the blob sizes in general) may be considered as 
adjacent so as to have iV= ( (4/3)ttL 3 | -1 . This factor disappears when substituting the mentioned volume for 
the £ integral. We thus arrive at: 



J"*-'— *..(?)H^i'> 



128tt 2 

\v\<L 



The last factor of the integrand points to an influence of the blob dimensions which is not accounted for 
by the autocorrelation function. It turns out that the effect of these dimensions L can only be neglected if L 
(being of the order of the scale of turbulence I) exceeds the maximal relevant value of \~?\. In view of the 
oscillating character of the other exponential factor this maximum is of the order of 



2Tr\ko\u sc — Ui\) =2tt ( 2k sin - J > 



& being the "scattering angle" between the propagation directions TQ of the incident radiation and Q P of 
the scattered energy. This angle being small in practice, we arrive at the condition Z>27r(/c tf) -1 , or #>X//, 
in order that the blob-size effect may be neglected. This inequality characterizes the region beyond that of 
the predominant forward scattering, that is the region of most interest for propagation up to distances well 
beyond the horizon of the transmitter. In this latter region it is also justified to extend the "u integration up 
to infinity. The conventional formula for the scattering coefficient, viz, 



(coVo) 2 (5e>2 



f 4e^oU ac -^)<7 5 Y|Y (12) 



128a- 2 

o<MO 

then results; it refers to scattering angles included by the unit vectors u sc and u%. 
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The relation between this expression and the power spectrum P(w) of the random fluctuations, to be defined 
with the aid of the three-dimensional Fourier transform of the spatial distribution of these fluctuations is well 
known. If P(«) is normalized such as to obtain the unit constant when integrated (over all real values of 
wi, a>2, and a> 3 ), we find 

We next consider the correlation, according to the Fraunhofer approximation, between the fields at two 
neighboring points Pi and P 2 . In view of the short distance of these points in practice, the y integral in (11) 
is hardly changed from its value for the previous case in which the angle \p vanished. However, the £ integration 
becomes more complicated, the vector V now being different from zero. Its evaluation indicates noticeable 
correlation between the fields at Pi and P 2 up to distances for which the angle \J/ = PiQoP 2 becomes of the order 
X//. In particular, the correlation is well established between points situated on one and the same line through 
the scattering source Qo, and least in directions perpendicular to such a line. 

6. Scattering Considered as a Markovian Process 

Statistical properties could in principle be derived for the higher-order scattering contribu- 
tions in a similar way as indicated for the Born approximation referring to first-order scatterings. 
The results would depend on the averages <C>fe(Qi) de(Q 2 ) . . . <5e($ n )> of multifold products. 
The corresponding final expressions holding for the combination of all terms of the Neumann 
expansion (2) would become extremely complicated though progress in this direction has been 
obtained [Hoffman, 1959, and Furutsu, 1963]. The difficulties here arising can be avoided 
when the successive scattering contributions, represented by the various terms of (2) , may be 
considered as independent of each other. The Markovian character of the scattering mecha- 
nism, then to be assumed, may be made plausible as follows. 

A narrow pencil of rays leaving some point source only produces a noticeable scattered 
field after it has covered some minimal distance. The order of magnitude of the latter can be 
estimated by investigating the scattered field produced by an antenna which radiates isotrop- 
ically, as given by (5), into an infinite random medium without boundaries. In the case of 
homogeneous turbulence the variance of this field follows rigorously from (6) by taking P 1 =P 2) 
the integrations extending over the entire space. Its computation has been worked out in 
[Fannin, 1956] assuming an isotropic auto-correlation function <<5e((?i) 8e(Q 2 )^> of gaussian 
form with respect to the distance Q1Q2', this function could also account for an anisotropy caused 
by wind effects. Fannin's main result (for zero wind velocity) can be interpreted thus that 
the variance of E {1) , divided by the squared primary field, does increase in proportion to d 3 / 
(d d cr 2 ) up to distances of the order of d cr , while tending to a limiting value of the order of d/do 
[see (8)] for d^>d cr ', the critical distance d cr =l 2 /\ is introduced here as a function of the scale 
of turbulence I and the wavelength X. The increase of the variance for small distances shows 
that the scattered field only becomes appreciable when the radiation has traveled over a distance 
of the order of dl /s d 2 J? if this quantity proves to be smaller than d cr , that is if d <Cd cr . This 
latter very general result proves to be independent of the special correlation-function model 
chosen by Fannin, as can be shown with the aid of (10). 

The importance of the dimensionless parameter d/d CT =d\/l 2 , here becoming obvious, can 
be understood as follows. Forward scattering is mainly restricted to a cone of angular devia- 
tions (around the undisturbed propagation direction) which are not larger than X/Z. An 
originally sharp ray having traveled over a distance d in a random medium is therefore broadened 
to a beam with a cross section having a dimension of the order of d\/l. This cross dimension 
just equals the size of a turbulent blob if d\/l=l, that is if d=d cr . Therefore, if d^>d cr , the 
radiation due to the scattering will have traversed, on the average, at least one blob. Also, 
d/dcr constitutes a measure of the number of blobs passed after traveling over a distance d in 
the forward direction. The geometric optical approximations of section 4 refer to small values 
of d/d cr . 
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In view of the above remarks individual scatterings take place, on the average, once along 
each section of a length of the order of rZJ /3 d 2 J r 3 if d <Cd C7 , but of d if d ^d cr . Hence the scattering 
angles associated with the successive individual scatterings may be considered as independent 
of each other if the correlation between the refractive-index fluctuations at the end of such a 
section proves to be negligible. This leads to the condition r/ 1/3 rf cr 2/3 >/ if d<^d cri and d o y>l if 
d^>d CT . It is therefore always sufficient to have d^lxi \<CJ (a condition satisfied for forward scatter- 
ing in order that the scattering angles connected with the nt\\ and (w-l)st scattering should 
be statistically independent. The angular spread caused by the nth scattering then only de- 
pends on the propagation direction obtained after the (n-l)st scattering, and not on those 
obtained after the preceding scatterings. This makes plausible the Markovian character of 
tropospheric multiple scattering (when occurring at all) . 

7. Equations Characterizing Markovian Small-Angle Scattering 

Physical phenomena associated with scattering generally depend on both angular and 
lateral deviations suffered by the energy bent aside by the irregularities of the medium. This 
suggests to introduce a joint probability density h z fixing the distributions of both types of 
deviations. Let us consider an energy unit (or a particle in the case of scattering of a stream 
of incident electrons or atomic particles by other particles having random positions in the 
medium traversed) which leaves an origin in the Z-direction along a path which would be 
rectilinear in the absence of scattering irregularities. We then define the function h z such that 
hz (X, ) r , x, y) dXdYdxdy constitutes the probability for the unit in question to pass the special 
plane Z=constant through a prescribed surface element dXdY while traveling there in a direc- 
tion comprised in an infinites timal cone, likewise prescribed. This cone subtends a solid-angle 
element dx dy if x and y are direction cosines with respect to the X and Y axes (both perpen- 
dicular to the Z-axis). The transport equation concerns the change of the function h z along a 
line element ds=(dX 2 +dY 2 +dZ 2 )*. 

Our analysis will be confined to forward scattering which implies that all relevant angular 
deviations (with respect to the Z-axis) may be considered as small. The limiting case of iso- 
tropic Rayleigh scattering is then excluded. The theory dealing with scattering angles of any 
magnitude usually resort at once to expansions in terms of Legendre functions, as discussed, 
e.g., [Lewis, 1950]. 

In the small-angle case the angle between two directions fixed by the cosines (xi, ?/]) and {x 2 , 2/2) can be 
approximated by 

i)~ sin#~{{ Xl -X2) 2J r(yi-y2) 2 }x. 

Moreover, all direction cosines with respect to the Z-axis may be replaced by unity. 

The small-angle approximations involve a considerable simplification. As an example of its significance 
we first consider the equation which expresses the property that scattering over a distance Z can be split into 
a pair of completely independent scatterings along any two subsections Z — Z\ and Z\, such in view of the 
Markovian hypothesis. The number of units passing the plane at a distance Z — Zi from 0, with lateral devia- 
tions £, 77, and angular deviations u, v, is determined by the quantity hz-z^, V, u, v). A fraction of these units 
will arrive, after traveling over a further forward distance Z\, within a prescribed range dXdYdxdy of lateral 
and angular deviations at the Z plane. During its propagation over Zi this fraction suffers a total lateral 
deviation 

\(X-Z l u-^+iY-Z.v-rj) 2 ]^ 

relative to the rectilinear path covered if it would have continued in the direction (u, v) up to the plane at the 
distance Z; in fact, along this latter path it would have arrived there at a point with coordinates Zi?/. + £, 
ZiV+y, Z. The corresponding angular deviation amounts to { (x — u) 2J r (y — v) 2 } 1 /}. In view of the axially 
symmetric character of the scattering, assumed throughout, the probability for these latter deviations would 
be the same as that in the case of lateral displacements X — Z\U— £ and Y—Z\V—t\ in the X and Y directions 
respectively, combined with angular displacements corresponding to a diviation away from the Z-axis into a 
direction with cosines x — u and y — v. Hence the probability for arriving within the range dXdYdxdy at the 
Z plane, after having passed the Z — Zi plane in the range d£ drj dudv, can be represented by 

hz^X — ZiU — %, Y—Ziv—tj, x—u, y — v) dX dY dx dy. 
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An integration over all possible situations in the (Z—Z{) plane, while taking account of the probability density 
hz-z^Z: V, u } v) for each of these situations, leads to the desired equation, viz, 

h z (X, Y,x,y) = \\\\ d£dvdudvhz-z 1 (£,v,u,v)hz 1 (X—£ — Z 1 u,Y--7i — ZiV,x—u, y — v). (13) 

— 00 

The integration over the direction cosines has been extended here up to infinity. This is allowed since the 
integrand will already be negligibly small near the limiting values ± 1 of the direction cosines with respect to 
the X and Y axis, which would correspond to very large angular deviations. Therefore, the integrations can 
just as well be continued over the complex directions fixed by cosines beyond these limiting values. Such 
infinite integration limits will be used throughout in our further analysis when the integration should comprise 
all real propagation directions. 

Any equation characterizing a special scattering medium should depend on its differential 
scattering cross section a(-y/x 2J ry 2 ), a dZ dx dy constituting the probability for first-order 
scattering, during propagation over a forward distance dZ, into an infinitesimal range of direc- 
tions (fixed by the cosines x and y) around a central axis in the Z direction. In view of the axial 
symmetry a only depends on the small angle (x 2J ry 2 ) H of this direction with the Z axis, but it 
proves to be convenient to define also a function <x(x } y) of two variables connected with that 
of one variable according to 

a(x, y)=(rQ x 2 +y 2 ). 

The final equation completely fixing all scattering phenomena can then be represented by 

h z (X, Y, x,y) 

(14) 



df = V dX +y W + dz) hz ^ Xy Y) x ' y ^ = | d * lv *& ^ hz ( X > Y > x -t, y-v)- 



The first two members of this special case of Boltzm arm's transport equation determine 
the change of h z along the line element ds, the coefficient in front of d/d z being unity in view of 
the small-angle assumption. The integral term represents the energy units that arrive at 
the plane Z= constant at prescribed points (X, F), traveling in directions prescribed by cosines 
x— £ and y — ??, and which are scattered near this plane into the direction fixed by the cosines 
x and y under consideration. The last term of (14) represents the loss due to the chance 
a(x, y)dxdyds for the units traveling along ds to be scattered aside into a new direction making 
an angle (x 2J ry 2 )* with the original direction. Extending the integrations up to infinity, the 
total amount of these losses becomes 

ds dxdy a(x,y)=-^> say. (15) 

Working out the integrations with the aid of (12), the quantity d defined here proves to 
be the same as that entering in the relation (8) for the varianceof the phase fluctuations. Accord- 
ing to (15) d represents the distance along which any energy unit becomes a chance 1 : 1 for 
being scattered aside completely. Therefore, this distance has been termed "mean free path 
for scattering." All terms of the transport equation (14) have thus been made clear. 

8. Solution of the Transport Equation in Terms of Laplace Transforms 

The integral term in the transport equation (14) constitutes a convolution product. 
This special form, which can only be obtained in the small-angle approximation, enables the 
solving of the equation with the aid of Laplace transforms; this has been worked out [Snyder 
and Scott, 1949] for the projection of the scattering paths on a plane through the Z-axis. In 
a completely four-dimensional treatment the transform in question may be defined by 
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Laiplh z (X } Y,x,y) = f ( f tdXdYdxdy h z (X, Y, x,y)e-*&-* T -**-**. 

Above we introduced the function h z by considering a unit element leaving the origin 
along the Z-axis. This implies the boundary condition h z = 5(X)5(Y)8(x)5(y) at Z=0. The 
Laplace transform of the solution satisfying this condition proves to be given (in so far as 
existing for special values of p h q h p 2 , #2) by: 

Lapl h z (X, Y, x, y) ==e' z/d o e d o J ° ' (16) 

as shown [Bremmer, 1963]. The function H(p,q) is defined here as the two-dimensional 
Laplace transform of <r(x,y). The function (16) also satisfies the Laplace transform of (13), as 
it should. 

The general solution obtained here contains as special cases the distribution functions 
connected with the lateral deviations only, or with the angular deviations only. The former, 
h z (X,Y) say, results after integrating h z (X,Y,x,y) over all possible values — °° <Cx,y<Cco of the 
direction cosines x and //. This integration corresponds to taking p 2 = q 2 =0 in the four-dimen- 
sional general Laplace transform. We thus arrive at the following Laplace transform for this 
function of two variables only: 

do I d d oJH(<hpJ, doqj) 

LaplA*(JT,r)=«-*/*« Jo . (17) 

The corresponding two-dimensional Laplace transform for the distribution function 
h z (x,y) of the angular deviations is obtained by substituting p 1 =q 1 =0 in (16). The final 
result reads: 

Laplhsfay^ezw**-™. (18) 

An integral expression equivalent to this relation has been derived in [Moliere, 1948]. 

9. Poisson Distribution Connected With Multiple Scattering 

The function h z includes the scattering contributions of any order. Each individual 
scattering being associated with the scattering coefficient a, the contribution of all scatterings 
of nth order can be recognized by the occurrence of n factors a. The role of this function is 
taken over by that of H in the Laplace transforms. Hence the expansion with respect to II 
of any result expressed in these transforms will show the distribution over the contributions 
connected with specific numbers of scatterings. 

As an example we consider all energy units that have passed through a plane Z = constant. The units 
applied in the definition of h z represent probabilities, hence their total number should equal unity. This 
total number is obtained, e.g., by integrating h z (x, y) over all values — °o <x, ?/< °° . The result equals the 
two-dimensional Laplace transform of a(x, y) at p = q = 0. Hence we should have 

1=1) dxdyh z (x, ?/) = {Lapl h z (x, y)} p = Q = Q , 

or, in view of (18), 

l = e z[mo, o)-Vd } =e -z/dQ y^ Z n \H(0, 0) \ n ^ 

The nth term then constitutes the fraction of units that has suffered n successive scatterings. Moreover, 

H(0,0) = l/do (19) 
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holds in view of the definition of d . The probability P(n) for n scatterings during forward propagation 
over a distance Z can therefore be represented by 

P( n )= e -z/d iKI^}l (20) 

This Poisson distribution here directly results as a consequence of the assumption of a Markovian scatter- 
ing mechanism. Equivalent derivations have been given by Dexter and Beeman [1949] and Fejer [1953]. 

10. Fokker-Planck Equation as an Approximation of the General Transport 

Equation 

The influence of the medium under consideration on the transport equation (14) results 
from the occurrence of the function <r(£, 77). However, the equation only proves to depend on 
the discrete set of parameters: 

rm idxdy a(x,y)(x 2 +y 2 ) n 

2ird dp <r(p) P *n + i = J-»J =Wj say? (21) 

I dxdya(x,y) 

which constitute the averages of the even powers of the scattering angles &i = (x 2J ry 2 ) 1/2 con- 
nected with a single scattering. In fact, a Taylor expansion of the integral term in (14), as 
worked out by Bremmer [1964], leads to the following alternative representation of the trans- 
port equation: 

(" ar+" w + £z) h ^ Y > x > *> =k S iw (S+UT *•<*■ Y > *> ti ■ ^ 

This equation suggests to consider the approximation obtained by restricting the right- 
hand side to its first term the significance of which has been discussed by Middleton [I960]. 
The approximation in question leads to a diffusion equation of the Fokker-Planck type, viz ; 



/ d 2 , d 2 \ , 4d Q ( d . d & \ . n 



Such an equation has been applied to particle scattering as early as 1929 [Bo the, 1929], at 
least insofar as the angular deviations are concerned. 

The exact solution satisfying this F. P. equation and also the previous boundary condition, that is h z 
= 8(X)8(Y)8(x)8(ij) for Z = } reads as follows: 

{h Z (X, y »3|y)lF.P. = y2( 02}2Z4 C W ' ( 23 ^ 

This expression is equivalent to a similar one derived [Rossi and Greisen, 1941] for scattering by particles. 
Obviously, its applicability requires sufficiently small values of the higher-order moments (#?), (#?), and so 
on. For the sake of completeness we also mention the corresponding probability densities for the lateral 
and angular deviations only: 

3rf 

We infer an increase of the average lateral deviations which is roughly in proportion with Z 3/2 , and of 
the angular deviations with Z l/2 . 
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1 1 . Evaluation of Scattered Fields for Special Antenna Systems 

Tropospheric radio propagation depends, among; other tilings, on the radiation diagrams 
of the transmitter and the receiver. That of the former may be given by the gain g tT U , *A) as 
a function of two independent angles fixing the directions of the rays leaving the transmitter. 
We shall define f as the angle between such a ray and the local horizontal plane through the 
transmitter, and \p as the angle between (a) the vertical plane through the ray in question, 
and (b) the other plane through the transmitter, the receiver and the center of the earth. 
The radiation diagram of the receiver can be fixed by the gain g r (£' , f) of the latter for radi- 
ations reaching it along directions fixed by two other angles f and \p' ; these angles may con- 
veniently be defined in an identical way with respect to the receiver. r I no analysis sketched 
below leads to the following ratio of the received and transmitted energies: 

^ j##fftrfo*)J^ (25) 



Ptr 16tt 

The lateral and angular deviations, away from the original direction fixed by f and \p< are completely 
fixed by the four angles f, \p, $*', and ,// when observing in a plane through the receiver thai is perpendicular 
to the original direction. We can determine the components X, Y, x, y of these deviations with respect to a 
coordinate system the X and Y axes of which are perpendicular to the (f, \p) ray leaving the transmitter, the 
A r axis being situated in the plane mentioned under (b). An elementary, though tedious geometrical analysis, 
shows that these components are given by 

respectively, a being the earth's radius and 8 the angular distance from the transmitter to the receiver. These 
latter expressions only result as approximations if #, f, \p, f ' and i// are assumed as small quantities of one and 
the same order of magnitude. The propagation distance covered in forward direction may then be taken equal 
to a iJ. 

The products d£d\p and df f d\p f constitute solid-angle (dements of infinitesimal pencils of rays leaving the 
transmitter or reaching the receiver respectively. Therefore, in view of (26), the probability for an energy 
unit to leave the transmitter along a given (f, yf) direction and to reach the receiver within a prescribed <l'<;'d\p f 
range is given by 

fc«*|a*(|+r)' atty, 0+r+f, tf+*'|ilr«(r% *')#'#', (27) 

in which 

^rec(f, V)=^~ (Jre C {{', +' ) 

constitutes the effective area of the receiver for the radiation arriving from the direction (£', \p'). 
The energy leaving the transmitter in a prescribed range d^dxp further amounts to 

The fraction of this energy that reaches the receiver within the other range d£'d\p' is obtained by multiplying 
the latter quantity by (27). The total received energy then results from an integration of this product over 
all relevant values of f, \p, f ' and \p', which leads to the final expression (25). 

The f, f r integrations may be restricted to positive values of f and f ; in fact, negative 
values correspond to rays suffering a reflection against the earth's surface while the effect of 
such rays could be included (at least for far-distance propagation) in the relevant radiation 
diagram. The conventional computations referring to the Born approximation would corre- 
spond to a substitution of the expression 



Ai ,} = 



*oa{x,y)KxY-yX)u(j) 



for h z , U(a) being Heaviside's unit function (unity or zero according as a>0 or a<0). 
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Apart from the difficulties connected with the exact determination of h z , the influence of 
the radiation diagrams will involve new mathematical complications. It is therefore recom- 
mendable to resort to simple representative situations. For instance we may consider the 
combination of a very narrow transmitting beam with an isotropic receiver. Under idealized 
conditions this corresponds to 

g r being independent of f ' and \p\ In this special situation the received field only depends on 
the lateral deviation a d 2 /2 away from the single original propagation direction £ = \f/ = 0. The 
resulting field is thus completely conditioned by the function h z (X,Y) of (17), and the general 
formula (25) here reduces to 

P rec X* (a$ 2 \ roQ . 

0*lV7 (28) 

An application of this relation is discussed in the next sections. 

12. Expansion for the Distribution of the Lateral Deviations Only and That of 

the Angular Deviations 

The last example shows the importance of the function h z (X, Y) in tropospheric far-distance 
propagation. According to (17) this function depends on the Laplace transform H(p, q) of 
the scattering coefficient <r(x, y). Assuming, once again, axially symmetric scattering, a merely 
depends on (x 2 +t/ 2 ) 1/2 and, consequently, H on (p 2 -\-q 2 ) 1/2 . For a very general class of scatter- 
ing models the latter dependence is given by an entire function. The corresponding Maclaurin 
expansion can be represented by 



H{p, g)= /(AW) 1 £ Cre (V^ +(f r, (29) 

a a /z =o 

in which c =l in view of (19). The introduction of the function f(s) of one variable enables a 
reduction of (17) to: 

ha,plh z (X,Y) = e- z/d oe <Wp2 ~ N2 . 

In view of the above mentioned Maclaurin series this relation proves to be equivalent 
to the following symbolical representation: 



LapU 2 (X,F) = e- z /"(N 






(30) 



Each power of b/dv here merely amounts to a multiplication by the same power of ZVp 2 +g 2 . 

The transform variables now only occur in the last exponential which, fortunately, can be 

recognized as the Laplace transform of an elementary function. In order to show this, we 

start from the relation: 

z,xv>+? 2 1 

-. Lapl . (31) 

-yjf+q 2 2iri-y]x 2 +y 2 -\ 2 

which holds when the following conditions are fulfilled simultaneously: 



(a) -yJx 2J ry 2 —\ 2 has to be defined with a positive real part, 

978 



(b) Vp 2 +# 2 has to be positive imaginary, 



(c) 3mX>0. 

This relation can be checked by introducing polar coordinates p, p in the double Laplace integral 
of its right-hand side. The <p integration then reduces to a zero-order Bessel function. The 
remaining p integral constitutes a special example of Sommerfeld's integral for the wave 
function of a point source; this integral represents the decomposition of this wave function 
into solutions of the wave equation that are separated in cylindrical coordinates. 

A differentiation of (31) with respect to X yields the relation to be applied to the last 
exponential in (30). We thus arrive at the following representation for the function h z itself: 



h z X,Y=- 



2*1 



,ife-} 



(X 2 +Y 2 -y 2 Z 2 f 



After a substitution of the Maclaurin series for J(s) we may work out an expansion of the 
exponential still occurring with respect to d/dy. We only need the odd powers, the even powers 
vanishing when applied to the subsequent function of 7. Each odd power of b/dy only affects 
a single term of the binomial expansion of this last function with respect to y 2 Z 2 /(X 2 -\-Y 2 ). 
An elementary, though tedious evaluation finally results in an expansion of the function h z 
itself which starts with: 



h z (X, Y) 



ZVc 



z 



.c A Z .c^Z* c; 

M rf 6 $ + 48 



27riL2^ (.Y 2 +F 2 ) 3/2 

I oorS c r>Z /CjC 4 . C 2 C 3 \Z 2 /CiCf CJC 3 \Z 3 

° \~6 <V V 10 ^ Ujdt + V.36 ^ 32 )d% 



8<%J 



ZJ 



(X 2 +Y 2 ) 512 



3 + elc 2 Z 



144 d% 3840 



lodlf(X 2 +Y 2 ) 7l2+ •' J' (32) 



In special applications, such as given below, the convergence of this series proves to be 
sufficiently rapid. 

The method, applied here, can be worked out in a similar way for the angular deviations, 
starting from (18) instead of (17). The series corresponding to (30) then proves to read: 



h z (x,y)-- 



Z 

d n 



•2wi\_{x 2 +y 2 y 



*4 



/ Z Z 2 J\Z* 

■" V' 3 d +ClCi d~l + Q d 



J) 

/ 



225 | c s 



Z ±( ^, , Z 2 , cM+CiC 3 ) Z 3 t clc 2 Z* 



c\ 



dl ' 6 dt ' 120 






(x 2 +y 2 ) 7 



13. Far-Distance Field in the Case of Norton's Model for Tropospheric 

Scattering 

The derivation of this field in the case of a narrow-beam transmitter and a wide-angle receiver is briefly 
indicated below. 

According to (28) the received energy is proportional to h z (X,Y) = h 7 (^ X 2 — Y 2 ) taken for the arguments 



= 4a 2 ' 



Substitution into the expansion (32) yields the following dependence on the linear propagation distance Z: 

+ 225 - 64 i3r7o + VX + ^J^ + Vl8+T6;^ + W^ + T^^|^ (33) 
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For c\9^ we get a dominating contribution, if Z> do, which decreases in proportion to Z -4 . This solution, 
corresponding to multiple scattering, involves a field strength decreasing as Z~ 2 . On the contrary, for ci = 0, 
the dominating term of h z decreases as Z -6 , and the field strength as Z~ z . A vanishing value of C\ occurs for 
the ''modified Bessel model," originally proposed by Muchmore [Norton and Norton, 1956], which is char- 
acterized by the isotropic autocorrelation function 

CUp1D = (p,1)Ki( p /1). 

This model is in accordance with many observational data [Norton, 1959], as well as the Villars and Weisskopf 
[1955] theory of turbulent mixing (turbulent blobs with dimensions conditioned by the vertical gradient of the 
refractive index). In the small-angle approximations, the modified Bessel model can be characterized by the 
isotropic case (/*= 1) of the derivation given in appendix 2 of [Norton, I960]. The scattering function is 
therefore given by 

-(«,„)- , 3 f W) ]M - (34) 

1 +7^x(^ 2 +^ 2 ) 



The corresponding function /(s) of (29) reads: 



^■(,-vf.) 



/(«) = « V 2 l-n/^r « • (35) 



This can be proved by putting (34) in the form: 

2B^ 2 d 2 1 



a(x, y) = 



irdo dB 2 (x 2 +y 2 +B)v 2 
for J3 = (tf?)/2, 



and by applying (31) thereafter in order to obtain the two-dimensional Laplace transform. 

The Maclaurin expansion of (35) yields the coefficients Ci. Their substitution into (33) results in 






For the transmitter and receiver here considered the decrease of the energy proves to be pro- 
portional to Z -6 in the plural-scattering range (r/>>c/ ) ; this corresponds to a field decrease 
proportional to Z -3 . 

14. Final Remarks 

The conventional theories for volume scattering in random continuous media only consider 
single scatterings (Born approximation). However, multiple-scattering effects may become 
relevant under special circumstances, as shown by Bugnolo [1960] and Anient [1960] for far- 
distance tropospheric radio propagation. It is therefore suggested to introduce in that case 
first of all a phenomenological theory taking account of the radiation diagrams of the transmitter 
and the receiver, without specifying the number of scatterings. Such a theory may start from a 
joint probability density for the lateral and angular deviations suffered by an energy unit 
which would travel along a rectilinear path in the absence of the random fluctuations of the 
medium properties. This probability density h z has to satisfy an integral equation, the "trans- 
port equation"; it reduces to a much simpler second-order differential equation of the Fokker- 
Planck type provided that very special approximations are justified. The general integral 
equation depends on the scattering coefficient introduced in the theories dealing with single 
scattering. The latter equation can be solved rigorously with the aid of Laplace transforms, 
the solution including the contributions associated with all orders of scattering. These con- 
tributions can be identified individually by a proper expansion of the solution. The Born 
approximation fits in this general theory by taking a proper function for h z . 

The application of the general theory to the modified Bessel scattering model, which agrees 
with the Villars and Weisskopf turbulence theory, leads to a decrease of the far-distance field 

980 



in the multiple-scattering range which is proportional to the inverse third power of the distance. 
This result is quite different from the exponential decrease conditioned by the Fokker-Planck 
approximation, or also From the corresponding decrease associated with exponential models 
(instead of the algebraic model of Norton) for the autocorrelation function fixing the scattering- 
function a. The above shows that the complete transport equation has to be applied in tropo- 
spheric propagation theories in order to get reliable results for the nonexponential autocorrela- 
tion functions confirmed by empirical data. 



The author acknowledges his contact with Dr. A. E. Goertz (CRPL Laboratory, Boulder, 
USA) who directed in particular his attention to the comprehensive literature on the subject. 
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